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Abstract 

We derive uniform approximations for contributions to Gutzwiller’s periodic-orbit sum for the 
spectral density which are valid close to bifurcations of periodic orbits in systems with mixed phase 
space. There, orbits lie close together and give collective contributions, while the individual contri¬ 
butions of Gutzwiller’s type would diverge at the bifurcation. New results for the tangent, the period 
doubling and the period tripling bifurcation are given. They are obtained by going beyond the local 
approximation and including higher order terms in the normal form of the action. The uniform ap¬ 
proximations obtained are tested on the kicked top and are found to be in excellent agreement with 
exact quantum results. 
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1 Introduction 


Semiclassical approximations in terms of periodic orbits belong to the main tools for the examination of 
spectral properties of quantum systems. They allow, for example, to explain fluctuations in quantum 
spectra in terms of the periodic orbits of the corresponding classical system. Semiclassical periodic- 
orbit approximations have been derived in cases where the classical dynamics is chaotic or integrable 
or has more general symmetries I, I, I, I 11, I, |. In these systems periodic orbits are typically 
either isolated or appear in families. 

Most systems, however, are neither chaotic nor integrable but show a complicated mixture of 
regular and chaotic behaviour. In these systems semiclassical approximations are more complicated 
and up to date there do not exist complete semiclassical approximations for quantities like the spectral 
density in terms of periodic orbits. These difficulties are due to the fact that in many situations 
periodic orbits neither appear in families nor can they be treated as being isolated in a semiclassical 
approximation. This is the case when there are other periodic orbits closely nearby. If the action 
differences of neighbouring periodic orbits is not large in comparison to h then the orbits yield a 
collective semiclassical contribution, and this is the typical situation when bifurcations of periodic 
orbits occur. In order to extend semiclassical approximations to systems with mixed phase space one 
has to derive expressions for the joint contribution of orbits which participate in a bifurcation. For 
generic two-dimensional systems this was first done by Ozorio de Almeida and Hannay. They derived 
an approximation which is valid in the vicinity of a bifurcation Q. In the present paper we extend 
the results of Ozorio de Almeida and Hannay and derive analytic formulas which interpolate over the 
regime from a bifurcation up to regions where the orbits can be considered isolated. 

Bifurcations are a characteristic phenomenon in systems with mixed phase space. They are respon¬ 
sible for the rapid increase of the number of periodic orbits when an integrable system is transformed 
into a chaotic system, e. g. by changing an external parameter. If one changes this parameter by an 
arbitrarily small but finite amount, then in general an infinite number of bifurcations occur, since 
they take place any time that the stability angle of a stable orbit is a rational multiple of 27r. There 
are different kinds of generic bifurcations, but the number of different forms is limited. The generic 
bifurcations that occur in two-dimensional conservative systems (or, equivalently, one-dimensional 
area-preserving maps) were classified by Meyer and Bruno |] 
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|. They are characterized by 
normal forms which describe the characteristic classical motion in the vicinity of a periodic orbit. 
Altogether one has to distinguish five different cases. These are the period-m-tupling bifurcations 
with m = 1 up to m = 5. They have the property that a central periodic orbit bifurcates and other 
periodic orbits split from the central orbit whose primitive period is m times the primitive period of 
the central orbit. (An exception is the case m = 1 for which there is no periodic orbit before the 
bifurcation.) The cases for m > 5 follow the same pattern as for m = 5. 

Ozorio de Almeida and Hannay derived their approximation by an integration over the normal 
forms. Their results are expressed in terms of diffraction catastrophe integrals. We call this approx¬ 
imation in the following the local approximation, since it is valid in the vicinity of a bifurcation. In 
farther distance from a bifurcation it reduces to an approximation of Gutzwiller’s type for isolated 
periodic orbits, it does not yield, however, the correct semiclassical amplitudes of the orbits. In |^] 
the results of Ozorio de Almeida were extended for the case m > 5 by including higher order terms 
in the normal forms. A slightly generalized version of this result is given in [M]. In the present paper 
we derive corrections for the low order bifurcations with m = 1 up to m = 3. We obtain uniform 
approximations in closed form which are valid from a bifurcation up to the regime where Gutzwiller’s 
approximation can be applied. 

The article is organized as follows. In section 2 we introduce the semiclassical method for treating 
bifurcations and in section 3 we present the results for the uniform approximations in two-dimensional 
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conservative systems that are derived in the appendices. In section 4 we give the corresponding results 
for one-dimensional area-preserving maps. Numerical examinations of the uniform approximations are 
carried out in section 5 on the example of a kicked top, and in section 6 we discuss the range of validity 
of our approximations and possible further extensions. 


2 Contributions to the spectral density 


In autonomous systems with discrete energy spectra the density of states 

d{E) = Y KE - En) = - -Im Tr G{E) 

^ TT 


( 1 ) 


can be expressed in terms of the trace of the (retarded) Green function. 

Semiclassical contributions of periodic orbits to the level density are derived from (||) by inserting 
the semiclassical approximation for the Green function. In a mixed coordinate-momentum represen¬ 
tation this approximation is given for two-dimensional systems by 


G{q',p,E) 


iTiy/2'Kih 


E 


- X '' X'K 


( 2 ) 


where the sum runs over all classical trajectories with initial momentum p and final position q' at 
energy E. The function for a trajectory ^ is defined as 


/ f^' 

S^{q,P,E)= p-dq + q-p, 

JQ 


( 3 ) 


where the integral is evaluated along the trajectory and q is the initial position of the trajectory 
which is determined by the initial momentum p and the final position q'. Its value follows from the 
condition that the gradient of the right-hand side of (^) with respect to q vanishes. S{q',p,E) is the 
generating function for a canonical transformation from final to initial coordinates of a trajectory. 
The mixed representation has been chosen since for a bifurcating orbit this transformation is close to 
the identity which cannot be generated in a pure coordinate- or momentum-representation. 

Finally, is the determinant of a matrix of second derivatives of 


( d'^St 


Dc = det 


dq'dp dq'dE 

d'^s^ a'^s^ 

\ OEdp EE^ 


( 4 ) 


and is an integer related to the Maslov index (see e. g. 
The semiclassical spectral density is then determined by 


d{E) = — ilm Jd'^q' d‘^pG{q',p,E) exp 


( 5 ) 


In the vicinity of a periodic orbit the integrals are evaluated by choosing local coordinates with 
one coordinate along the orbit and one perpendicular to it. If the integral over the perpendicular 
direction is evaluated in stationary phase approximation one obtains the Gutzwiller approximation 
for the semiclassical contribution of the orbit. 
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Here is the action of the orbit whose value is given by the function S^{q',p, E) — q' ■ p for the 
considered periodic orbit. T^, and are the period, stability matrix and Maslov index of the 

orbit, respectively. The integer denotes the repetition number of the orbit. In this notation we 
consider a periodic orbit and its multiple traversals which all give semiclassical contributions to the 
level density as different periodic orbits. The uniform approximations for bifurcating periodic orbits 
that are derived in this paper will be expressed in terms of exactly the same classical quantities that 
appear in (^). 

We consider now the contribution of a bifurcating orbit with repetition number r to the level 
density. The condition for such a bifurcation is that the stability matrix of the corresponding primitive 
periodic orbit (repetition number r = 1) satisfies TrM = 2cos(27m/r) with integer r and n and thus 
the stability matrix of the r-th traversal has a trace which is equal to two. Let I be the greatest 
common divisor of r and n. Then the bifurcation is a period-m-tupling bifurcation with m = r/l. 
Near the bifurcation the integral over the perpendicular coordinates cannot be evaluated in stationary 
phase approximation since there are stationary points nearby which correspond to the other orbits 
participating in the bifurcation. This is reflected by the fact that the Gutzwiller approximation 
(^) diverges at the bifurcation. Instead one has to derive a joint contribution of all orbits which are 
involved in the bifurcation. This is achieved by expanding the generating function S{q',p, E) in higher 
order around the central orbit. In general this results in a complicated exponent in the integrand of 
(^). The integrals can be considerable simplified by a canonical transformation of the coordinates and 
by using the fact that the form of eq. dH) is semiclassically invariant under canonical transformations. 
This follows from work of Miller |15] and is discussed by Littlejohn |16]. 

The most simple form that the generating function can take near the bifurcation is given by the 
normal form. This normal form contains the information about the number and arrangement of the 
orbits which are involved in the bifurcation. The transformation to the normal form coordinates 
has the further advantage that then the integral over the coordinate along the periodic orbit can be 


performed trivially. These steps are explained in detail in [13| and for that reason we give here only 
the resulting formula 


d^{E) 


1 




Re 


dq' 


^ 1 dS 

d'^s 


dpdq' 


eyiV{^S{q\p,E)-^(^p-^-v] . (7) 


Here p and q' are coordinates in the Poincare surface of section perpendicular to the orbit, and 
S{q',p, E) is the generating function for the r-th iterate of the Poincare map which obeys the conditions 


dS , dS dS ^ 


( 8 ) 


where T is the time from initial to final point. 

The approximation of Ozorio de Almeida and Hannay for the contributions of orbits near a bifur¬ 
cation is obtained by inserting the normal form of the generating function S{q',p,E) for a particular 
generic bifurcation. This yields the semiclassical contributions in terms of standard diffraction catas¬ 
trophe integrals. For the bifurcations which are considered in this paper the normal forms are given 
by 

m = 1 : S{q',p,E) = So{E) + q'p — ^p^ — eq' — aq'^ 

m = 2: S{q',p,E) = So{E) + q'p-^p"^ - eq'^ - aq"^ (g) 

m = 3; S{q',p,E) = SoiE) + q'p - ^{q''^ + p"^) - ^{p^ - 3pq'^) , 

where e is a parameter which is zero at the bifurcation and cr is a sign factor. As will be shown in 
the appendices, the diffraction integrals for these cases can be expressed in terms of Bessel functions. 
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The properties of the bifurcations corresponding to these normal forms are discussed in more detail 
in the next section. All these bifurcations involve two periodic orbits. 

The approximation of Ozorio de Almeida and Hannay is valid in the vicinity of a bifurcation and 
hence is called local approximation by us. In farther distance from a bifurcation the local approxima¬ 
tion splits up asymptotically into a sum of separate contributions of Gutzwiller’s type for the periodic 
orbits. However, in this limit the semiclassical amplitudes of the orbits have a fixed ratio, i. e. there 
is a fixed relationship between the stabilities (and periods) of the different orbits. In more detail, the 
approximation holds when the following relations between the monodromy matrices and the periods 
of the orbits are valid. 


m = 1 : Tr Ml Tr M 2 - 4 = 0 , Ti = T 2 

m = 2 : Tr Ml-|-2 Tr Mo — 6 = 0 , Ti = Tq (10) 

m = 3 : Tr Mi -|- 3 Tr Mq — 8 = 0, Ti = Tq . 

For m = 2 and m = 3 the index 0 denotes the central bifurcating orbit and the index 1 the satellite 
orbit. For the case m = 1 there is no central periodic orbit and the two periodic orbits which are 
involved in the bifurcation are given the indices 1 and 2. The relations (0) follow from the normal 
forms and are valid in the vicinity of a bifurcation. In farther distance from a bifurcation, however, 
they are not valid any more, and then the local approximation becomes inaccurate. In order to obtain 
a formula which uniformly interpolates over the region from the bifurcation up to regimes where 
the Gutzwiller approximation is valid (without restrictions on the semiclassical amplitudes) one has 
to include higher order corrections to the normal forms in (P). The resulting integrals can then be 
reduced to simpler forms by appropriate coordinate transformations. These calculations are carried 
out in appendices Si and P, and the results are discussed in the next section. 


3 Results for the uniform approximations 


3.1 The isochronuous bifurcation 

The bifurcation which is described by the normal form in (^) for m = 1 has the following property: 
On one side of the bifurcation where e and a have opposite signs there exist two periodic orbits, one 
stable and one unstable. We denote these orbits by and ^ 2 - On the other side of the bifurcation 
where e and a have the same sign both orbits are complex, and we give the index 1 to that orbit for 
which the imaginary part of the action is positive. Due to the shape of the function S{q',p, E) — q'p 
this bifurcation is called tangent bifurcation or saddle-node bifurcation. 

The uniform approximation for the semiclassical contributions of the two orbits is derived in 
appendix Its form is different on the two sides of the bifurcation. On the side where the orbits are 
real it is given by 


d^{E) = 


1 

nh 


27rA5 

sn 


—sign(AS') 







V h 


( 11 ) 


.-.3(^)-..3(^V 


which is invariant under exchange of the two indices. The quantities which appear in (^) are the 
mean action S = (^i -|- S2)/2 , the action difference A5 = (^i — «S'2)/2 and the mean Maslov index 
= { 1^1 + ^ 2 )!“^ = 12 + cj/2 of the orbits. Here v is the index in (|^ and a is the sign factor in the 
normal form (^. Furthermore, Aj denotes here and in the following the semiclassical amplitude of an 
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orbit 


( 12 ) 


* M, - 2) ’ 

where rj = sign (Re (TrMj — 2)). Properties of the two orbits when they are real are listed in table 1, 
and an expansion of the classical properties of the orbits in terms of the coefficients in the (extended) 
normal form are given in appendix 


m = 1 

a > 0 

(7 < 0 

ae < 0 

unstable, Uu = 

stable. Vs = V + 1 

Ss>Su 

unstable, = v 

^s stable, Vg = V — 1 

Su>Ss 


Table 1: Properties of real orbits that participate in a generic isochronuous bifurcation (m = 1). 
The orbits are denoted by and where the indices u and s denote the unstable and stable orbit, 
respectively. 

On the other side of the bifurcation where the orbits are complex the semiclassical contribution to 
the level density is given by 


d^{E) = Re 
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nh 


2AS 


ttH 


1/2 


I - ITT ITT 

—b - v - -a 

h 2 A 


Ai + Ao 
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exp 

f\AS\\ , A,-A2 




+ 




(13) 


where now the actions and amplitudes are complex. Again is the index in (^) and a is the sign 
factor in the normal form i). Both equations (jl^) and (H) can be written in a combined form by 
expressing them in terms of an Airy function. 

In the limit e —> 0 the leading order semiclassical contribution of equations (pT|) and (IB) is given 
by 

, ToFfl/S) / So TT TT \ , 

d^(E) = cos[^ - -V - -a] . (14) 


aAl'^ V ^ 2 2 ) 

Here Tq and Sq are, respectively, the period and action of the orbits at the bifurcation, a is the 
coefficient in the normal form (0), and I is the repetition number of the bifurcating orbits. The 
contribution (14) is of order larger than the contribution of an isolated period orbit. All 

appearing classical quantities in (^) depend on the integer 1. In detail, Tq^/ = lTo^i=i, Sq^i = lSo^i=i, 
ui = and ai = The relation for the Maslov index follows from the Maslov index of 

the unstable orbit in table 1 for multiple traversals. The last relation is obtained by noting that all 
coefficient in the normal (^ increase by a factor I as I is increased. However, a is restricted to be ±1. 
This is achieved by a subsequent (canonical) transformation p ^ p/'/i and q —> qVl. It follows that 
the amplitude of the contribution at the bifurcation decreases like with 1. 


3.2 The period-doubling bifurcation 

The period-doubling bifurcation which is described by the normal form in (|^ for m = 2 has the 
following form: On one side of the bifurcation where e and a have the same sign there is only one 
orbit which is called the central orbit. This orbit changes its stability at the bifurcation from stable 
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to unstable or vice versa, and a new orbit appears which is named the satellite orbit since the two 
fixed points of the Poincare map which belong to this orbit lie symmetrically on both sides of the fixed 
point of the central orbit. This bifurcation is also called pitchfork bifurcation. 

The uniform approximation for the contributions of these orbits to the level density is derived in 
appendix ^ and is given by 


d^{E) = Re 


1 

vrh 

Ai 


ttAS 


2% 

+ — 

2 ^v/2. 


1/2 


I - nx iTX 
exp I — D - —V - rcr 


(15) 


+ 




This approximation is valid on both sides of the bifurcation. Here S = {Si + So)/2, AS = {Si — So)/2, 
V is the index in (^) and a is the sign factor in the normal form (^). The values of v and a can be 
determined from the properties of the orbits which are listed in table 2. Furthermore cJi = sign(AS') 
and 02 is a sign factor which discriminates between both sides of the bifurcation. It is 1 when both 
orbits are real and —1 when only the central orbit is real. The orbit contributes also to (15) when 
it is complex, but its action and amplitude factor are always real. The expression of the classical 
properties of the two orbits in terms of the coefficients in the (extended) normal form are given in 
appendix ^ 


m = 2 

a > 0 

cr < 0 

a > 0 

e > 0 : ^0 stable, vq = + 1 

e < 0 : ^0 unstable, vo = ^ 

stable, 1^1 = -|- 1 

Si>So 

e > 0 : ^0 unstable, 1^0 = i' 

e < 0 : ^0 stable, i 2 q = u — 1 

unstable, 121 = 1 / 
Si>So 

a < 0 

e < 0 : .^0 unstable, vq = 

e > 0 : .^0 stable, vq = v + 1 

unstable, 1^1 = 
Si<So 

e < 0 : ^0 stable, i/q = 1 / — 1 

e > 0 : ^0 unstable, 1/0 = 1 / 

stable, 1/1 = 1 / — 1 

Si<So 


Table 2: Properties of real orbits that participate in a generic period-doubling bifurcation (m = 2). 
The central orbit is denoted by and the satellite orbit by .^i. 


In the limit e —> 0 the leading order semiclassical contribution of equation (15) is given by 


Tor(l/4) 


47r/v/^fi5/4^i/4 [ n 


So vr 


TT 


7T 




(16) 


Here Tq and Sq are, respectively, the period and action of the orbits at the bifurcation, a is the 
coefficient in the normal form (|^, and I is the repetition number of the orbit The contribution 
(16) is of order larger than the contribution of an isolated period orbit. The dependence of the 


classical quantities in (IT) on the integer I are given by: Tq,; = ITq^i=i, Sq^i = lSo,i=i, vi = 11 ^ 1=1 and 
ai = l^ai=i. The relation for the Maslov index follows from the Maslov index of the unstable orbit 
in table 2 for multiple traversals, and the relation for the coefficient a from considerations analogous 
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to those for the isochronuous bifurcation. It follows that the amplitude of the contribution at the 
bifurcation decreases like with 1. 


3.3 The period-tripling bifurcation 


The period-tripling bifurcation is described by the normal form in (^) for m = 3. It involves two orbits, 
the central orbit .^o and the satellite orbit Both orbits exist before and after the bifurcation. As 
e goes through 0 both orbits approach each other, they coincide at the bifurcation and then separate 
again. For that reason this bifurcation has also been named “touch and go” bifurcation |^]. 

The uniform approximation for this bifurcation is derived in appendix and is given by 


, I /27r|A5| f z - iTT 

= -ReJJ^exp -S--. 


^0 \ 
2 +2^/3; 
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J_i/6 


^0 

2 




/|AS| 

f\^S\\ 




(17) 


Here S = (Si + So)/2, AS = {Si — So)/2, ly = i^i is the Maslov index of the satellite orbit. Furthermore 
a = sign(A5). Properties of the two orbits are listed in table 3, and the expression of the classical 
properties of the two orbits in terms of the coefficients in the (extended) normal form are given in 
appendix 


m = 3 

£ > 0 

e < 0 


^0 stable, uq = u + 1 

unstable, 1 / 1 = 1 / 

Si<So 

^0 stable, i/Q = 1 / — 1 

unstable, 1/1 = 1 / 

Si>So 


Table 3; Properties of orbits that participate in a generic period-tripling bifurcation (m = 3). The 
central orbit is denoted by and the satellite orbit by . 


In the limit e —> 0 the leading order semiclassical contribution of equation ( p^ is given by 


, /pN _ ^or(i/6) 

^ 9/7r3/2h"/3a2/3 


cos 


TT 

n ~ 2 ^ 


(18) 


Here Tq and Sq are, respectively, the period and action of the orbits at the bifurcation, a is the 
coefficient in the normal form (|^, and I is the repetition number of the orbit ^i. The contribution 
is of order larger than the contribution of an isolated period orbit. The dependence of the 


classical quantities in (1^) on the integer I are given by: Tq^i = ITq^i=i, Sq^i = /S'o,/=i, t'/ = li^i=i and 
a; = lai=i. The relation for the Maslov index follows from the Maslov index of the unstable orbit in 
table 3 for multiple traversals. It follows that the amplitude of the contribution at the bifurcation 
decreases like Z“^/3 1. 


4 Uniform approximations for maps 

Though we concentrated on autonomous systems in the preceding paragraphs, the uniform contribu¬ 
tions given there can equally be applied, with minor modifications, to quantum maps. These maps 
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are described by a unitary time-evolution operator F, the so-called Floquet operator; the dynamics 
of the system is generated by repeated applications of the operator on a state in Hilbert space. F 
has eigenstates and unimodular eigenvalues where the phases (pi are called quasi energies and 

the states quasi-stationary states since many quantum maps originate from a stroboscopic description 
of periodically driven systems. For such systems with Hamiltonian H{t + T) = F[{t), the unitary 
time-evolution operator after period T is chosen as F = U{T), implying U{nT) = F^. 

Restricting oneself onto a finite-dimensional subspace of Hilbert space with dimension N the quasi¬ 
energy spectrum can be obtained by solving the secular equation 

N 

P(x = e-*^) := det(F - A/) = ^ = 0 

n=0 


for the N x A^-matrix F. Here the set of traces TV F'^ comes into play: their knowledge for n up to N/2 
allows to construct the first half of the coefficients a„ via the so-called Newton formulae; the other 


half follows from the unitarity of F which entails self-inversiveness of the secular polynomial |15]. 
The semiclassical starting point for the calculation of TV F’^ is an expression which is nearly identical 
its appearance with equation for autonomous systems, 


m 


TtF*" = 


dpdq 

2Trh 


d^S 


dpdq' 


exp 


\^^{S{q\p\nT) - q'p) - 


(19) 


The major difference to the expression for autonomous systems is that one does not switch to energy 
representation by a Fourier transform with respect to time t. In translating the contributions which 
were derived for autonomous systems into the appropriate contributions to TV F’"' attention has to be 
paid to the following five difference: (i) The orbits which contribute are those with a fixed period n, 
not those with a given energy E] (ii) the primitive periods have to be expressed in units of T and thus 
are integer valued; (hi) the action is not the reduced energy dependent one, but depends on time (that 
is, on the number n); (iv) instead of taking twice the real part, the full complex contribution has to 
be taken; (v) the results further differ by a factor 2'Kti. Since these are mainly formal differences the 
morphology of the contributions remains unaltered. 

In the limit —> 0 each orbit of primitive period no contributes individually according to 


TVF” 


E 

orbits 


^0 

2-TVM|V2 


exp 




( 20 ) 


In the following we denote the contributions to TVF” by (7^”^ and use Ai = no,i(^?(TrMj — 2))“^/2 
abbreviate the stability amplitudes with p = sign (Re (TV Mi — 2)). S = {Si Si)/2 is the mean action 
and AS = {Si — Si)/2 the action difference where z is 0 or 2 depending on the considered bifurcation, 
and u is the index in (pT|). 

For m = 1 the collective contribution is 
1/2 


C(n) ^ 


2-kAS 


3h 


Ai A2 


.S 


.vr 


exp 


^- 1/3 




+ J] 


1/3 


(^) 


( 21 ) 


—sign(A5) exp 




when both orbits are real. When they are complex one obtains 




2AS 


Trh 


1/2 


exp 


l ^ ITT ^ 7 ^ 

— S - v - -a 
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The amplitudes Ai are now complex quantities, a is the sign factor in the normal form A pair of 
orbits involved in a period doubling bifurcation [m = 2) gives the contribution 


Cin) ^ 


ttAS 

2h 


1/2 

exp 



iaiir/S 




(23) 


with (Ti = sign (A5) and cj 2 = 1 when the satellite is real, (T 2 = — 1 otherwise. For m = 3 the result 
reads 


cin) ^ 


' 27 r|AS'| f / - zvr 


^0 \ 
2 +2^/3y 




(24) 


with a = sign (AS). 

These expressions give also the contribution of repetitions of the bifurcating orbits if one substitutes 
the corresponding classical quantities. For the l-th repetition, the stability angle w of a stable orbit 
in TrM = 2 cos a; and the instability exponent u of an unstable orbit in TrM = ±2 cosh u increase 
linearly, uoi = loJi=i and ui = lui=i, as does the action Si = lSi=i and the Maslov index of the unstable 
orbits The Maslov index of the stable orbits is + sign(5;^^^ — 5^“^). 


5 Numerical results 


In this section we want to test the quality of the uniform collective contributions that we derived for 
the various types of bifurcations on a certain quantum map, a periodically kicked top |jT^, ^]. It 


will turn out that the uniform approximations indeed are accurate both close to the bifurcation as 
well as in far distance; the local approximation is only valid close to the bifurcation while far away 
the orbits can be treated as being isolated via the stationary-phase approximation. 

Tops are dynamical systems that involve the angular momentum operators Jx, Jy, Jz, satisfying 
the usual commutation relations [Jk, Ji] = i^kimJm, where T% is set to unity. The evolution of the 
system is such that the total angular momentum + Jy + J'^ = j{j + 1) is conserved. This introduces 
the well known good quantum number j which fixes the Hilbert space dimension 2j -|- 1. j -|- 1/2 
further plays the role of the inverse of Planck’s constant; accordingly, the semiclassical limit is reached 
by sending j ^ oo. 

The specific top used here is described by the Floquet operator 


F = exp(-/ 2 ^.^^ 1 - ipzJz) exp{-ipyjy) 


X Jx) 


(25) 


The dynamics consists of rotations by angles pi and nonlinear rotations (torsions) of strength ki. For 
the study of bifurcations we hold the pi at fixed values {px = 0.3, py = 1.0, pz = 0.8) and take k = 
as our control parameter, with kx = fe/10. The classical counterpart of the system is integrable at 
/c = 0 and displays well developed chaos at A: = 5. 
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There is a very convenient testing tool which enables one to examine each contribution of a given 
cluster of periodic orbits individually. It is obtained by considering the function 

1 Jmax 

= - ^e“*^'^TrF^(j) 

Jmax jmin +1 

J —Jmin 

In its essence this is a Fourier coefficient of Tr F” with respect to j; finite limits have to be taken for 
practical reasons if one wants to evaluate the sum for the quantum system, they further give control 
over the desired degree of rigor of the semiclassical limit. We used jmin = 1 and jmax = 64 with a 
single exception for m = 3. 


(26) 
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Figure 1: Sequence of action spectra |r(^)(5)p as the control parameter crosses a tangent bifurcation 
{m = 1) at fe = 2.45. A ghost peak shows up already before the bifurcation at S' = 5.3; beyond the 
bifurcation the peak splits into two, each corresponding to a bifurcating orbit. 


peak height 




Figure 2: Height of the peaks in the action spectrum at the positions of the two orbits that are 
engaged in the bifurcation of Figure |^. Both the uniform and the local bifurcation work well close to 
and on the left of the bifurcation. The local approximation starts to fail as k is increased, and the sum 
of isolated contributions of Gutzwiller’s type becomes valid. Note that the uniform approximation 
remains accurate and can hardly be distinguished from the exact result in this plot. 


From the asymptotic behaviour (2C) of all uniform approximations it is clear that one expects 
peaks in at values of the argument which correspond to actions of periodic orbits. The 

function is called the action spectrum for that reason. Figure || confirms that peaks indeed 
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show up. Plotted is a sequence of action spectra for n = 1 as the control parameter k is varied across 
a tangent bifurcation at k = 2.45. Already before the pair of orbits comes into existence a peak is 
visible at 5 = 5.3. This peak arises from the complex predecessors of the bifurcating orbits that were 
already observed in |^|. The other peaks pertain to different orbits. Slightly beyond the bifurcation 
the new-born orbits have nearly identical actions and give rise to a single peak that would be resolved 
if one would go to much higher jmax- Increasing the control parameter further the peak splits into two 
peaks which are located at the now well separated values of the orbits’ actions. 

The quality of our results can now be tested by calculating the height of the peaks |T("')(5 ci)P both 
quantum-mechanically exact as well as on grounds of the various semiclassical approximations. Figure 
^ depicts the height of the aforementioned peaks as the control parameter is varied across the tangent 
bifurcation. In the vicinity of the bifurcation both the local as well as the uniform approximation are 
accurate to a degree that makes it difficult to resolve the error at all. The uniform approximation 
remains also valid as one moves away from the bifurcation. There the sum of individual contributions 
of isolated orbits is seen to gain validity. 


action spectrum 



Figure 3: The exact peak in the action spectrum of Figure [T| at fe = 3.0 is compared to semiclassically 
evaluated ones. The uniform approximation works well. The local approximation predicts peaks of 
equal height. Treating both orbits as isolated gives peaks that are slightly too high. 

Figure ^ illustrates the superiority of the uniform approximation over the local one away from 
the bifurcation. To that end, the exact peak itself is compared to the approximated ones. The local 
approximation assumes that the stabilities of both orbits are equal and gives two peaks of same height. 
With the sum of two isolated terms of Gutzwiller’s type both peak heights are slightly overestimated. 
The uniform approximation almost coincides with the exact result. 

In Figure |^, the peak heights of orbits participating in a period doubling at A: = 4.3 are plotted. 
In this case it is hardly possible to discriminate between the local and the uniform approximations. 
Before the amplitudes start to differ significantly, erratic deviations to the exact result are encountered 
which stem from overlapping peaks of other orbits. One such orbit has already been included here 
with its isolated contribution for k > 5.8 and cures part of the problem, but other orbits become 
relevant from k ~ 6.2. Instead of further improving on this we settle for the information that the 
approximations work well close to the bifurcation and turn to other period doublings. Figure ^ actually 
arises from a similar situation of two overlapping peaks; both peaks, however, arise from two pairs of 
orbits which undergo period doublings at A; = 2.78 and k = 2.84, respectively. That both pairs have 
almost the same action is a pure coincidence. Plotted is now the collective contribution of the total 
number of four orbits in various semiclassical approximations and the exact result. The situation is 
qualitatively the same as for the tangent bifurcation: both the local and the uniform approximation 
are excellent close to the bifurcations; moving away the uniform approximation remains valid and the 
sum of isolated orbits starts to work well while the local approximation breaks down. 
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Figure 4: Peak height as in Figure for orbits involved in a period doubling bifurcation (m = 2) at 
k = 4.3. The local and uniform approximations are again accurate before and close to the bifurcation. 
This time they remain so until overlapping peaks from other orbits give rise to deviations. One such 
orbit has already been included, others not, as is explained in the text. The sum of individual terms 
that treats all orbits as isolated fails close to the bifurcation and regains validity far away. 


peak height 




Figure 5: Two pairs of orbits involved in period doublings at k = 2.78 and k = 2.84 give rise to 
overlapping peaks in the action spectrum; the height at the actions of two bifurcating orbits is plotted 
here. The accuracy of the uniform approximation in the whole parameter range, the failure of the 
isolated treatment close the bifurcation, and the breakdown of the local approximation at a distance 
to the bifurcation is clearly visible. 


peak height 



Figure 6: Height of peaks in the action spectrum that arise from orbits involved in a period tripling 
bifurcation (m = 3) at A: = 3.545. The satellite is engaged in a tangent bifurcation with another 
satellite at A; = 3.525; this requires to increase jmin = 2^^ + 1 as discussed in the text. The additional 
satellite is only felt slightly at the bifurcation and the accuracy is high for the local and uniform 
approximation. 
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Figure 7: The same period tripling is investigated as in Figure For the chosen = 1 an additional 
satellite overshadows the result close to the tripling. As one goes away from the bifurcation to higher 
k the uniform and isolated approximations work fine and the local approximation is inaccurate. 


For the period tripling we hit on difficulties from the close neighborhood of a tangent bifurcation 
of the satellite. This tangent bifurcation was found to be very close to the tripling for all investigated 
cases. The problem can be overcome if jmin is chosen large enough in order to ensure that the further 
satellite which is involved into the tangent bifurcation is well separated, the separation being measured 
by its difference in action in units of Planck’s constant, (j + 1/2)A5. Then our analytical expression 
(^), which does not describe the tangent bifurcation, can be tested without being overshadowed by 
the existence of the additional satellite. Since the relevant action difference is tiny we raised jmin up 
to 2^^ + 1, with jmax = 2^^ + 64. The result is presented in Figure For the lower values of k the 
tangent bifurcation at A; = 3.525 is still felt strongly and all semiclassical approximations go wrong. 
At the period tripling (k = 3.545) and beyond it the uniform and the local approximation work well 
again. The breakdown of the local approximation is observed if one goes to a different scale of the 
control parameter, as shown in Figure |^. 

In conclusion, for all studied cases the uniform approximation was found to give excellent results (as 
long as the results were not overshadowed by the existence of further orbits). The local approximation 
is of the same quality close to the bifurcation. It breaks down when the relations (0) between the 
monodromy matrices become invalid as one increasingly detunes the control parameter. Along the 
same line, the action difference of the orbits in units of Planck’s constant is large in that region of 
parameter space, and the use of the sum over isolated orbits (^) makes sense there. 


6 Conclusions 

In this paper we examined semiclassical contributions of three different bifurcations in two-dimensional 
conservative systems and one-dimensional area-preserving maps. These bifurcations are the generic 
period-m-tupling bifurcations with m = 1, 2 and 3. We extended the local approximation of Ozorio 
de Almeida and Hannay for these cases and derived uniform approximations which interpolate over 
the regime from the bifurcation up to regions where the orbits can be considered isolated. The 
approximations were tested numerically on the example of a kicked top and were found to be in 
excellent agreement with exact quantum calculations. The local approximation was found to be good 
near a bifurcation, but its quality decreased in farther distance from a bifurcation. In contrast, 
the Gutzwiller approximation is good sufficiently far away from a bifurcation, but it diverges at a 
bifurcation. 

We derived the classical properties of periodic orbits that follow from normal forms with higher 
order corrections. We further examined the semiclassical contributions of the orbits directly at the 
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bifurcation (e = 0). It was found that with increasing m bifurcations contribute more strongly to the 
level density. In more detail, the semiclassical amplitude of the orbits at the bifurcation is of order 
for m = 1, 0{h~^/^) for m = 2, and for m = 3. Furthermore, for larger m the 

semiclassical amplitude decreases more slowly for higher repetitions of the bifurcating periodic orbits. 
For the l-th. multiple of the bifurcating orbits the semiclassical amplitude decreases by a factor 
for m = 1, for m = 2, and for m = 3. For the case m > 5 it follows from |13| that the 
semiclassical amplitude at the bifurcation is of order and it decreases by a factor with 

increasing 1. Thus bifurcations with larger m have a stronger influence on the level density. 

All uniform approximations in this paper are valid as long as the participating periodic orbits do 
not bifurcate further, i. e. as long as the bifurcations can be considered isolated. It can happen that 
a periodic orbit undergoes several subsequent bifurcations. This has been observed in the numerical 
examples where a tangent bifurcation occurred very close to a period-tripling bifurcation. Similar 
situations for period-doubling and period-quadrupling bifurcations have been observed in H, 01- 
In these cases the uniform approximation has to be modified by fully integrating over an extended 
normal form in which higher order terms are included that describe also subsequent bifurcations. For 
the mentioned period-tripling bifurcation this is done in |^|. 

An extension of the present results concerns the case m = 4, i. e. the period-quadrupling bifurcation. 
Since the case m > 4 was treated in []^| this is the only remaining generic bifurcation. The treatment 
of the generic period-quadrupling bifurcation is more complicated than the other cases since it involves 
three periodic orbits whose action differences all scale with the same leading power in e (for small e). 
As a consequence one has to add two correction terms to the local approximation in order to obtain a 
formula which has the correct Gutzwiller limit. Furthermore, the diffraction catastrophe integral for 
this case can in general not be expressed in terms of simple functions. A treatment of this bifurcation 
is in preparation p^ . 

A further extension concerns systems with symmetries. In these systems there can be further 
kinds of bifurcations which are characteristic for the considered symmetry. Normal forms for sys¬ 
tems with symmetries have been derived in |^, ^|. For some cases the uniform approximations for 
bifurcations in systems with symmetries can be obtained from a slight modification of the formulas 
for generic bifurcations. For example, in systems with time-reversal and reflection symmetries there 
can be isochronuous pitchfork bifurcations. These bifurcations are described by formula (^) with 
the replacement of Ai by 2Ai. Another example where the generic formula for m > 5 is applied to 
symmetric period-n-tupling bifurcations with n = m/2 is discussed in [^. 
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A Uniform approximation for the isochronuous bifurcation 

The starting point for the derivation of the uniform approximation is the expansion of the generating 
function in the vicinity of the bifurcation 

S(q',p, E) = So + q'p - eq' - aq’^ - bq‘^ - . (27) 

Here we went one order higher in q' than in the normal form expansion, and in the following we will 
treat the extra term as a perturbation to obtain the uniform approximation. By rescaling q' and p we 
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can enforce a = ± 1 . The fixed points of the map generated by S are given by Sqi = p, Sp = q', and 
up to second order in \/[e[ we get for the two fixed points engaged in the bifurcation 


Pi,2 = 0 


(28) 


while we will not consider the additional fixed point which now formally arises as a third solution of 
the equations. The value of the actions S = S — q'p at the fixed points is 


5i,2 = 5o T 


2e 


£ 

3a 




+ 0(|e|5/2), 


(29) 


where Sq is the value of the generating function at the origin. If e has the same sign as a then the 
orbits are complex and there is an ambiguity for choosing the sign of the square root in Eqs. (|2^ ) and 
(^^. We then choose to give the index 1 to that orbit for which the imaginary part of the action is 
positive. This formally corresponds to choosing \/—ej (3a) = —isign(e)A/|e/(3a)|. 

The periods of the orbits are given by 


7"i,2 — 



+ £*(£) , 


(30) 


where Tq = dS/dE, evaluated at the origin, and corresponds to the mean period. The traces of the 
monodromy matrices follow from the relation 


™ = (^) 


-1 


1 + 


d'^S d^S 525525 ' 


dp dq' dp dq' dp^ dq'‘^ 
which is evaluated at the stationary points and leads to 

Tr Ml ,2 = 2 =F + 0{\£\^/‘^) . 

Furthermore the semiclassical amplitudes are given by 

7"i,2 _ 1 


^ 1,2 = 


l^?7i,2(TrMi,2 -2) ^|12ae| 


1/4 


ToT 


2b ^ \ I r 

^To + se) 


0 ( 8 ) 


(31) 


(32) 


(33) 


where pi = sign (Re (TrMj — 2)), and / = ri = r 2 is the repetition number of the orbits. 
In the following we will use the definitions 


5 = 


5i + 52 


A5 = 


S 1 -S 2 


A = 


+ T 2 


AT = 


Ti — Ay 


E = 


El + E2 


(34) 


where ei and E 2 are the Maslov indices of the orbits when they are real. 

We continue now with the evaluation of the integral in (^. The main contribution to the integral 
over q' comes from the region near the stationary points. For that reason we consider q' in the following 
as a quantity of order 0{\£\^/‘^). Then the exponent of the integral can be simplified by substituting 


.eh 96 

q = —^ + X — X — 
^ 9a2 3a 


which leads to a reduction of the generating function 


S{q',p,E) — q'p = S — £x — ax^ — —p"^ + 0{\£\^^'^) . 


(35) 


(36) 
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This is just the usual normal form with the replacement of Sq by S. Furthermore the exponential 
prefactor in the integral is modified by the Jacobian of the transformation. After an integration over 
p we obtain 


d^{E) = Re 


( i Q ITT,. in 

exp 

lV2^ 


dx (Tq — — eEx \ exp (^^x + ax^'^ 


, (37) 


where the exponential prefactor has been expanded up to 0{x). The integral in can be split into 
two terms. The first term with the constant term in the exponential prefactor has the same form form 
as the local approximation, and the second term is proportional to the derivative of the first term 
with respect to e. The integrals can be found in the section on Airy functions in |^|. The result 
depends on whether £ has the same or the opposite sign as a. We first consider the case when they 
have opposite sign, i. e. when the orbits are real. The contribution to the level density is then given 
by 


d^{E) = Re- 


exp 


( in i'E ^, in ^ A / _ 

2 4^) j 27rTo 




V 


e 

3a 


J. 




l-nie 

9|a| 


+ 


2hE 


3a 


J_ 


-'/'V n 

f\^s\\ 


2/3 


h 




(38) 


This result can be expressed also by the Airy function Ai and its derivative, but it seems more natural 
to express it by Bessel functions since then the analogy to the results for higher repetition numbers 
is more visible. All coefficients can be expressed by the classical actions and amplitudes. The final 
result is 




27rA5 


3h 


/ 3 / 5* TT 

.4COS 




—sign(AS') AA cos 


f\^S\\ 


(39) 



-T 


2/3 




Here we replaced v + a{2 by the mean Maslov index v. That this is correct can be seen from a 
stationary phase evaluation of the integrals which shows the vi ^2 = e + a/2 ± sign(A5)/2. The first 


term in the curly brackets in (^) is just the local approximation in which the stability factors of both 
orbits are equal. At the bifurcation (e —> 0) the dependence of the classical quantities on e guarantees 
a finite result. The second term in the curly brackets is the uniform correction which ensures the 
correct limit as h ^ 0 for finite e, since then one arrives at a sum of two individual contributions for 
the two orbits, each of Gutzwiller’s type. 

On the other side of the bifurcation where the orbits are complex and e and a have the same sign, 
the evaluation of the integral in leads to 


d^{E) = Re- 


exp 


in in ,, in ^ 


lV27r^h^ 





2i£ 


h J 3\/3|( 


£e + 


2f>To 

3a 




(40) 

In this equation the mean quantities S and A are real, but the differences AS and AA are purely 
imaginary. The quantities are defined in Eqs. (p^), and (p^) with the convention for the sign of 
the square root that has been discussed above 


dAE) = Re — 
vrn 


2AS 


eTi 


1/2 


exp 


l ^ IE IE 

—b - v - -a 

h 2 A 




(41) 
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In the limit that the argument of the if-Bessel functions is large one obtains only the contribution of 
that orbit for which the imaginary part of the action is positive 

d^(E) = Re ^ exp (^is, - . (42) 

This is connected to Stokes’ phenomenon. 


B Uniform approximation for the period-doubling bifurcation 


For the period-doubling bifurcation, we start again with the normal form S{q',p, E) of the action and 
go beyond the local approximation by incorporating the next order term in the expansion with respect 
to q' 

S{q',p, E)=So + q'p - eq'^ - aq'^ - bq'^ - ^p^ , (43) 

where once more u = ±1. No term of order q'^ can be present since both satellite hxed points belong 
to the same orbit and thus must have the same action. The fixed points lie at 


Po = 0, 
Pi = 0, 


% = o 


(44) 


q'l = ±\ - 


2a 




where the subscript 0 indicates the central orbit which has a repetition number of 21, and the subscript 
1 denotes the satellite orbit. The actions of the orbits are Sq and 



(45) 

dS/dE evaluated at the origin and 


Ti = To + + 0(e2) . 

2a 

(46) 


The traces of the monodromy matrix follow from (31) and are given by 

Tr Mq = 2 — 2as 

Q L 

Tr Ml = 2 + Aae - + O(e^) , 

and the stability prefactors follow as 

x4o = 


(47) 


x4i - \Aef\ V2(^:ro + ge + ^e + 0(e2)) , 


(48) 


where the repetition numbers of the orbits are vq = 21 and ri = 1. 
In the following we will use the definitions 


5 = 


Si + So 


AS = 


Si -Sn 


u = 


+ Vo 


(49) 


2 ’ 2 ’ 2 

In order to evaluate the integral (|^) we want to get rid of the dependence of the exponent on g'® 
and introduce for this purpose a new variable by 


/2 2 4 

g =x -X — . 

2a 


( 50 ) 
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If one considers q' as a quantity of order 0(|e|^/^) the generating function reduces to 


S — q'p = Sq — ex^ — , 


( 51 ) 


where a = a — be/{2a). One can check that the new prefactor in front of x'^ yields the correct action 
(^) up to order e^ at the fixed points. After integrating over p one arrives at 


d^{E) = Re 


exp f (j) 

/(27r7i)3/2 


J frr ^bTo 2 2^ f * / 2 I ~ 4^ 

ax I io-^ “ ^ex I exp [ex + ax 1 


(52) 

where the exponential prefactor has been expanded up to order 0{x‘^). Once more the integral splits 
into two terms of type 


= 


J dx exp — (ex^ — ax^)^ 


(53) 


e 

2d 


1/2 


exp -AS J_i /4 


A5 




AS 


Jaiw/S 


and 

B2 = 


J dxx^ exp -^(ex^ + ax^)^ 


(54) 


= -—Bi - — 

4a 8 a 


e 

2d 


1/2 


exp <{ -AS \ ( J 3/4 


AS 


^i(Tl7r3/8 


— 0'i5'2 J_3/4 


A5 


iCTi7r3/8 


where AS as given above, ui = sign(a) = sign(A5), and o '2 = sign(e). Now we introduce (T 2 = —crid '2 
to discriminate between both sides of the bifurcation (cj 2 = 1 when the satellite orbit is real and 
CJ 2 = — 1 when it is complex). Expressing all coefficients by the actions So,i and stability factors Agq 
the final result is 


d^{E) = Re 


1 

nh 

Ai 


ttAS 


2h 

^0 


1/2 


exp 


i - m ITT 

—S - V - -a 

h 2 A 


(55) 


2 


'2 V2 


CT2 Ji/4 (^) + J_i/4 




In contrast to the tangent bifurcation, all appearing classical quantities are real even when the satellite 
orbit is complex. In the limit that the argument of the Bessel functions is large the expression reduces 
to a sum of the two Gutzwiller contributions of the orbits when cj 2 = 1, and to the single Gutzwiller 
contribution of the central orbit when CJ 2 = — 1 . 


C Uniform approximation for the period-tripling bifurcation 

For the period-tripling bifurcation the calculations are more involved, since the orbits which participate 
in the bifurcation lie now in a plane. For the previous two bifurcations (m = 1 and m = 2) the orbits 
lay on a line so that one had to treat only one-dimensional integrals. 

The normal form of the generating function S{q',p,E) in terms of q'- and p-coordinates is given 
by 

S{p, q', E) = pq - + q'"^) - - 3m'^) • (56) 
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The function S{q',p,E) — q'p has four stationary points, one at the origin which corresponds to the 
central orbit, and three others which lie on an equilateral triangle and correspond to the satellite 
orbit. In contrast to the previous two cases, one has to go two orders beyond the normal form in 
the expansion of the generating function S since the first correction does not change the amplitude 
prefactor of the satellite orbit. The higher order terms in this expansion are restricted by the condition 
that the three stationary points which correspond to the satellite orbit have to yield all three the same 
action and the same trace of the monodromy matrix of the satellite orbit. This leads to the expansion 


S{p,q',E) = 


{p^ - 3pq'^) - ^(/ + + q'^) - ^p{Sp‘^q' - q'^) 


^ {p^ - 2p^q'‘^ - Spq'"^) - - 29pq'^) - ^^{-Sp'^q' +p‘^q'^) .(57) 


4v/2 


64v/2 


4\/2 


The requirement that the three stationary points yield the same classical properties allows also a 
further term of the form dq'{p^ + but the properties of the orbits do not depend on this 

term and it can be removed in a later step by modifying the coefficients (^) of the transformation 
(32). For that reason we do not write it. 

In the following the expansions of several quantities in terms of e are given. The action of the 
satellite orbit is 


Si = So- 


27a? 


e^- 


166_^4 _ 4(64b^ - 243a^ -24ac) ^5 ^ 


81a4 


729a6 


(58) 


and its period is of the form 


rr.rr.dS 


= n- 


Aee ,2 




9a2 


e" + 0{£^) . 


(59) 


The traces of the stability matrices of the periodic orbits follow from (^) and are given by 
TrMo = 2 — £^ 

rr.rr 2 (2511a'^ - 1286^ + I92ac) 4 

TVMi = 2 + 36^ + ^--- + 0{£^) . 

54a^ 


(60) 


The amplitude prefactors follow as 
To To 


= 


Ai = 


d/v'ITrMo - 2| 3/|e| 

Ti _ 1 

VITrMi -2| \/3/|e| 


Tn- 


Aee 


e — 


ro(2511a^ - 12852 + 192oc) 
324a4 


£2 + 0(e2) , (61) 


where the repetition numbers of the orbits are tq = 2>l and ri = 1. 

We continue now with the evaluation of the integrals in (j^ with the generating function of (57). 
The exponent in the integral is simplified by applying a transformation of the form 


p 


= p + cip^ + C2pq + cag^ + c^p^ + c^p^q + cepg^ + cyg^ 


g = g + dip + d 2 pq + d^q^ + d^p'^ + d^p q + d^pq + djq'^ 


(62) 


which removes the terms proportional to q''^p‘^ and q'^p^ ^ in the exponent (within the considered 
order of the approximation). The coefficients of the transformation are given by 


Cl = 


bV2 V2 


6a 


216a3 


(285^ - 243a^ - 24ac)£ 


C2 = 


3V2a 


C3 = -Cl 
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( 63 ) 


27a^ - 252 + 2ac 
~ 12^2 ’ C5- 


4 ’ 


C6 = 


351a^ - 1662 + I6ac 


di = ^ , d 2 = -2ci , ds = 0 , ^4 = -6 , 

27a^ - 1662 + I 60 C ,b 

- --' *-2’ 


96a2 


81a^ — 862 _|_ 
48a2 


C7 = 


'4 ’ 


and the transformation leads to 
d^{E) = 


2 rcxD 

Re 


/ oo roo 

dp / dq[ai + a2{p^+ f)] 

-OO J —OO 




(64) 


where the exponential prefactor in (|^ times the Jacobian of the transformation (^) has been expanded 
up to order e'^. The new constants appearing in 


are 


«! — To, 


0(2 


144£i;a2 + 251iroa^ - 128ro62 + 192roac _ 27a‘^l / Ap 


288a2 


4|£| V2\/3 


26 243a^ - 2862 + 24ac o 

d — a — —£ H-——o- £ 

3a 54a3 


(65) 


After a change of variables p = cos 4* and q = sin 4> the integrals in (|^ ) are transformed into 

1 pOO ^ 

~ sy J d^[ai + 202 /] exp (^Sq — el — a/^/2 cos ( 34 .)) I 

dl 


d^{E) = 


rRe 


Jo 

1 


Slirh"^ 

where the relation 


Re 


POO 

/ dl [oi + 202 /] Jo 
Jo 




p27r 


(i4> exp{iz cos(m4>)} = 27rJo(2:) 


( 66 ) 

(67) 


has been used. The two remaining integrals are evaluated according to appendix ^ and result in 


/o V ^ 
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f\^S\\ 
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+ iaJ, 


5/6 
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,(69) 


where A5 = (^i — So)/2 = —2£^/{27dJ) and a = sign(A5) = —sign(e). Altogether one obtains 
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where S = {Si + So)/2. 


D The diffraction integral for the period-tripling bifurcation 


In this section we evaluate the diffraction integral which appears in the uniform approximation for the 
period-tripling bifurcation. We consider first the case 2 ; > 0: 
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(73) 

(74) 
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The parameter e has been introduced in ( [7l| ) in order to make the integrals absolutely convergent. 
The integrals leading to (H) after a substitution x = can be found in |2^. From ([7^ to ([7^ the 
limit e —> 0 has been performed and the duplication formula of the Gamma function has been used. 
From ([73|) to (^) the sum has been split into three parts (by taking every third term, respectively) 


where the third part vanished. Furthermore the recurrence formula and the triplication formula of 
the Gamma function have been used. From (|7^ ) to (^) the definition of Kummer’s function iFi has 
been used, and from (^) to (^) the formula exp(z 2 )Ji,( 2 ;) = ( 2 / 2 )^ iFi(i/ -|- 1/2; 2v -|- 1; 2iz)lT{v -|- 1) 

[!]• 

The corresponding expression for negative values of 2 ; follows from the evenness and oddness of 
the real and imaginary part of the integral, respectively. We obtain 
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and from the derivative of this integral with respect to 2 ; follows 
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